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In the context of neutron star mergers, we study the gravitational wave spectrum of the merger 
remnant using numerical relativity simulations. Postmerger spectra are characterized by a main 
peak frequency /2 related to the particular structure and dynamics of the remnant hot hypermassive 
neutron star. We show that /2 is correlated with the tidal coupling constant K 2 that characterizes the 
binary tidal interactions during the late-inspiral-merger. The relation 42 ( 1 ^ 2 ) depends very weakly 
on the binary total mass, mass-ratio, equation of state, and thermal effects. This observation opens 
up the possibility of developing a model of the gravitational spectrum of every merger unifying the 
late-inspiral and postmerger descriptions. 

PACS numbers: 04.25.D-, 04.30.Db, 95.30.Sf, 95.30.Lz, 97.60.Jd 98.62.Mw 


Introduction.- Direct gravitational wave (GW) ob¬ 
servations of binary neutron stars (BNS) late-inspiral, 
merger and postmerger by ground-based GW intefer- 
ometric experiments can lead to the strongest con¬ 
straints on the equation of state (EOS) of matter at 
supranuclear densities m- There are two ways to set 
such constraint^ (I) measure the binary phase during 
the last minutes of coalescence using matched filtered 
searches [J |3H5] ; (11) measure the postmerger GW spec¬ 
trum frequencies using burst searches [iiil. 

Method (I) relies on the availability of waveform mod¬ 
els that include tidal effects and are accurate up to 
merger usiini- Here, “up to merger” indicates the end 
of chirping signal in a precise sense that will be described 
below. Tidal interactions are significant during the late 
stages of coalescence at GW frequencies /gw ^ 400 Hz 
(for typical binary masses), and affect the phase evolu¬ 
tion of the binary. The zero-temperature EOS is con¬ 
strained by the measure of the quadrupolar tidal cou¬ 
pling constant (or equivalent/correlated parameters, 
e.g. m that accounts for the magnitude of the tidal in¬ 
teractions mm- 

Combining results from numerical relativity and the 
effective-one-body (EOB) approach to the general rela¬ 
tivistic two-body problem HMS], one can show that the 
merger dynamics of every irrotational binary is charac¬ 
terized by the value of [H] . At sufficiently small sep¬ 

arations, the relevant dependency of the dimensionless 
GW frequency on the EOS, binary mass, and mass-ratio 
is completely encoded in the tidal coupling constanl|^ A 
tidal effective-one-body model compatible with numeri¬ 
cal relativity data up to merger was introduced in 
but no prescription is available to extend the model to 
the postmerger. 


^ GW observations of BNS mergers can also constrain the source 
redshift Eng. 

^ The spin dependence is approximately linear for small spins 
aligned with the orbital angular momentum. 


Method (H) relies on the high-frequency GW spec¬ 
trum, and can, in principle, deliver a measure indepen¬ 
dent on (I) [7]. Binary configurations with total mass 
M < Mthr ~ 2.9 Mq are expected to produce a merger 
remnant composed of a hot massive/hypermassive neu¬ 
tron star. The merger remnant has a characteristic GW 
spectrum composed of a few broad peaks around /gw 
1.8 — 4 kHz. The key observation here is that the main 
peak frequencies of the postmerger spectrum strongly 
correlate with properties (radius at a hducial mass, com¬ 
pactness, etc.) of a zero-temperature spherical equilib¬ 
rium star in an EOS-independent way [SI dH]. Thus, 
a measure of the peak frequency constraints the corre¬ 
lated star parameter. Recently, there has been intense 
research on this topic, and various EOS-independent re¬ 
lations were proposed lidHHIll. Most of the relations 
are constructed for equal-mass configurations and do not 
describe generic configurations for different total masses 
and mass-ratios, e.g. [H US]. Additionally, the post¬ 
merger GW spectrum might be influenced in a compli¬ 
cated way by thermal effects, magnetohydrodynamical 
instabilities and dissipative processes. 

In this paper we observe that the coupling constant 
can also be used to determine the main features of the 
postmerger GW spectrum in an EOS-independent way 
and for generic binary configurations, notably also in the 
unequal-mass case. The observation opens up the possi¬ 
bility of modeling the complete GW spectrum of neutron 
star mergers unifying the late-inspiral and postmerger 
descriptions. Geometrical units c = G = 1 are employed 
throughout this article, unless otherwise stated. We use 
/ for the spectrum frequencies and w for the instanta¬ 
neous, time-dependent frequency. 

Numerical Relativity GW Spectra .— The numerical 
relativity data used in this work were previously com¬ 
puted in [IlllSI- In our simulations we solve Einstein 
equations using the Z4c formulation m and general rel¬ 
ativistic hydrodynamics [28j . Our numerical methods are 
detailed in [2^ . The binary configurations consid- 
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FIG. 1. Simulations of BNS and GWs. Top: real part and amplitude of the GW mode Rh 22 /{vM) and the associated 
dimensionless frequency MtU 22 versus the mass-normalized retarded time t/M for a hducial configuration, H4-135135. The 
signal is shifted to the moment of merger, fmrg, defined by the amplitude’s peak (end of chirping). Also shown is (twice) the 
dynamical frequency MQ = dEb/dj ~ MljJ22 /2. Bottom: Snapshots of logj^Q p on the orbital plane, during the late inspiral 
(left), at simulation time corresponding to tmrg (middle), during the postmerger (right). 


ered here are listed in Table [T] In the following we sum¬ 
marize the main features of the GW radiation obtained 
by BNS simulations. 

We consider equal and unequal masses configurations, 
different total masses, and a large variation of zero- 
temperature EOSs parametrized by piecewise polytropic 
fits [31]. Thermal effects are simulated with an additive 
thermal contribution in the pressure in a T-law form, 
^th = (Eth — l)pe, where Tth = 1.75, p is the rest-mass 
density and e the specific internal energy of the fluid, 
see ISIIMIIIS]. The initial configurations are prepared 
in quasicircular orbits assuming the fluid is irrotational. 

Initial data are evolved for several orbits, during 
merger and in the postmerger phase for > 30 millisec¬ 
onds. A detailed discussion of the merger properties de¬ 
termined by different EOSs, mass, and mass-ratio is pre¬ 
sented in mus]. The binary configurations in our sam¬ 
ple do not promptly collapse to a black hole after merger, 
but form either a stable massive neutron star (MNS) or 
an unstable hypermassive neutron star (HMNS), which 
collapses on a dynamical timescale row « 


200 ms |37j. Both HMNS and MNS remnants at for¬ 
mation are hot, differentially rotating, nonaxisymmetric, 
highly dynamical two-cores structures, e.g. [33113H] . 

The typical GW signal computed in our simulations 
is shown in Fig. [^ for a fiducial configuration. We 
plot the real part and amplitude of the dominant I = 
m = 2 multipole of the s = —2 spin-weighted spheri¬ 
cal harmonics decomposition of the GW, R{h+ — ihx) = 
m — 2Yim{S,4)), versus the retarded time, t. The 
figure’s main panel also shows the £ = m = 2 instanta¬ 
neous and dimensionless GW frequency Mu >22 = Md(j)/dt 
where (j) = —arg(i?/i 22 )- The bottom panels show snap¬ 
shots of logj^Q p on the orbital plane, corresponding to 
three representative simulations times. 

The waveform at early times is characterized by the 
well-known chirping signal; frequency and amplitude 
monotonically increase in time. The GW frequency 
reaches typical values wgw = Stt/qw ~ ^22 ^ 0.1 — 
0.2/M, i.e. /gw ~ 0.8 - 1.6 kHz for a M = 2.7Mq bi¬ 
nary. The chirping signal ends at the amplitude peak, 
max \ Rh 22 \, which is marked in the figure by the middle 






































































































































3 



/ [kHz] 

FIG. 2. GWs spectra from BNS. The plot shows only a rep¬ 
resentative subset of the configurations of Table [I] Triangles 
mark frequencies /mrg corresponding to tmrg, circles mark /2 
frequencies. 


vertical line. We formally define this time as the moment 
of merger, t„irg, and refer to the signal at t > tmrg as the 
postmerger signal. The GW postmerger signal is essen¬ 
tially generated by the m = 2 structure of the remnant, 
see bottom right panel of Fig. The frequency increases 
monotonically to Mw 22 0.2—0.5 as the HMNS becomes 
more compact and eventually approaches the collapse. 
Assuming the remnant can be instantaneously approxi¬ 
mated by a perturbed differentially rotating star [38] , the 
/-mode of pulsation is strongly excited at formation and 
it is the most efficient emission channel for GWs. 

The GW spectra are shown in Fig. for a representa¬ 
tive subset of configurations. Triangles mark frequencies 
/mrg corresponding to tmrg- Circles mark the main post¬ 
merger peak frequencies /2 ~ 1.8 — 4 kHz. The small fre¬ 
quency cut-off is artificial and related to the small binary 
separation of the initial data; physical spectra monoton¬ 
ically extend to lower frequencies. From the figure one 
also observes that: (i) there exists other peaks, expected 
by nonlinear mode coupling or other hydrodynamical in¬ 
teractions [23l ISSl |38] ; (ii) peaks are broad, reflecting the 
nontrivial time-evolution of the frequencies (see Fig. 
and also the spectrogram in [35]); (hi) secondary peaks 
are present in most of the configurations, their physical 
interpretation has been discussed in [22l [231 |26l |38| . We 
postpone the analysis of these features to future work. 
In the following we focus only on the /2 peak, which is 
the most robust and understood feature of the GW post¬ 
merger spectrum. 

Characterization of the postmerger GW spectra .— 
Here, we show that /2 correlates with the tidal coupling 


TABLE I. BNS configurations and data. Golumns: name, 
EOS, binary total mass M, mass ratio q, /2 frequency in 
kHz, dimensionless M /2 frequency, tidal coupling constant 
K 2 ■ Gonfigurations marked with * are stable MNS. 


Name 

EOS 

M [Mo] 

q 

/2 [kHz] 

M/2 [xlO^] 

1. II 

SLy-135135 

SLy 

2.70 

1.00 

3.48 

4.628 

74 

SLy-145125 

SLy 

2.70 

1.16 

3.42 

4.548 

75 

ENG-135135 

ENG 

2.70 

1.00 

2.86 

3.803 

91 

SLy-140120 

SLy 

2.60 

1.17 

3.05 

3.906 

96 

MPAl-135135 

MPAl 

2.70 

1.00 

2.57 

3.418 

115 

SLy-140110 

SLy 

2.50 

1.27 

2.79 

3.426 

126 

ALF2-135135 

ALF2 

2.70 

1.00 

2.73 

3.630 

138 

ALF2-145125 

ALF2 

2.70 

1.16 

2.66 

3.537 

140 

H4-135135 

H4 

2.70 

1.00 

2.50 

3.325 

211 

H4-145125 

H4 

2.70 

1.16 

2.36 

3.138 

212 

ALF2-140110 

ALF2 

2.50 

1.27 

2.38 

2.931 

216 

MSlb-135135* 

MSlb 

2.70 

1.00 

2.00 

2.660 

290 

MSl-135135* 

MSI 

2.70 

1.00 

1.95 

2.593 

327 

MSl-145125* 

MSI 

2.70 

1.16 

2.06 

2.740 

331 

2H-135135* 

2H 

2.70 

1.00 

1.87 

2.561 

439 

MSlb-140110* 

MSlb 

2.50 

1.27 

2.08 

2.487 

441 

MSlb-150100* 

MSlb 

2.50 

1.50 

1.87 

2.303 

461 


TABLE II. Fit coefficients of different quantities at tmrg and 
of M /2 with the template in <§• 


Q(«r) 

Qo 

ni [xlO^] 

n 2 [xlO®] 

di [xlO^] 

^mrg 

-0.1201 

+2.9905 

-1.3665 

+6.7484 

j-mrg 

+2.8077 

+4.0302 

+0.7538 

+3.1956 


+0.3596 

+2.4384 

-1.7167 

+6.8865 

M/2 

+0.053850 

+0.087434 

0 

+0.45500 


constant that parametrizes the binary tidal inter¬ 
actions and waveforms during the late-inspiral-merger. 
The relation / 2 (tt|") depends very weakly on the binary 
total mass, mass-ratio, and EOS. We use a large data 
sample of 99 points including the data of mM- 

Let us first briefly summarize the definition of K 2 and 
its role in the merger dynamics. 

Within the EOB framework, tidal interactions are 
described by an additive correction A^lr) to the ra¬ 
dial, Schwarzschild-like metric potential A{r) of the EOB 
Hamiltonian m- The potential A{r) represents the bi¬ 
nary interaction energy. In order to understand its phys¬ 
ical meaning, it is sufficient to consider the Newtonian 
limit of the EOB Hamiltonian, Heob ^ Mc^ + + 

^{A{r) - 1) + 0(c-2), where ^ = MaMbUMa + Mb) 
is the binary reduced mass, p the momenta, and A(r) = 
1 — - = 1 — 2 -§M— with rAB the relative distance be- 
tween the stars (constants c and G are re-introduced for 
clarity). The tidal correction {r) is parametrized by 
a multipolar set of relativistic tidal coupling constants 
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FIG. 3. M /2 dimensionless frequency as a function of the tidal coupling constant . Each panel shows the same dataset; the 
color code in each panel indicates the different values of binary mass (top left), EOS (top right), mass-ratio (bottom left), and 
Fth (bottom right). The black solid line is our fit (see Eq. ([^ and Table]^; the grey area marks the 95% confidence interval. 


, where A, B label the stars in the binary mm- 
The leading-order contribution to {r) is proportional 
to the quadrupolar (/ = 2) coupling constants, nf = 
[Xa/CaY’ Mb/Ma where Ma is the mass of star A, 
Ca the compactness, Xa = Ma/M, and the i = 2 
dimensionless Love number [401115] . The total £ = 2 cou¬ 
pling constant is defined as -|- , and can be 

written as 


K 


T 

2 


= 2 


\{i+qrc% 


{i + qrci,j 


( 1 ) 


assuming q = Ma/Mb > 1. The leading-order term of 
the tidal potential is simply A^(r) = — 

A consequence of the latter expression for A^(r) is 
that the merger dynamics is essentially determined by 
the value of kJ [TO] . All the dynamical quantities develop 
a nontrivial dependence on as the binary interaction 
becomes tidally dominated. The characterization of the 
merger dynamics via is “universal” in the sense that 
it does not require any other parameter such as EOS, M, 
and q. (There is, however, a dependency on the stars 
spins.) For example, at the reference point tmrg) the cor¬ 
responding binary reduced binding energy the re¬ 
duced angular momentum and the GW frequency 


Muj^ 2 ^ can be fitted to simple rational polynomials [TO] 


Q(k^) = Qo 


1 -I- niK2 + 

1 -t- d\ ^2 


( 2 ) 


with fit coefficients {ni,di) given in Table [TT] 

In view of these results, it appears natural to investi¬ 
gate the depedency of the postmerger spectrum on 
Our main result is summarized in Fig. [^ which shows 
the postmerger main peak dimensionless frequency M /2 
as a function of kJ for a very large sample of bina¬ 
ries. Together with our data we include those tabu¬ 
lated in dUM]. The complete dataset spans the ranges 
M G [ 2 . 45 M 0 , 2 . 9 M 0 ], q G [1.0,1.5], and a large varia¬ 
tion of EOSs. The peak location is typically determined 
within an accuracy of Sf ^ ±0.2 kHz, see also [18]. Each 
of the four panels of Fig.|^shows the same data; the color 
code in each panel indicates different values of M (top 
left), EOS (top right), q (bottom left), and Tth (bottom 
right). The data correlate rather well with . As indi¬ 
cated by the colors and different panels, the scattering of 
the data does not correlate with variations of M, EOS, 9 , 
Tth- The black solid line is our best fit to Eq. where 
we set 712 = 0 and fit also for Qq, see Table [ll[ The fit 
95% confidence interval is shown as a gray shaded area 
in Fig. § 
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We argue that the observed postmerger correlation 
with is a direct consequence of the merger univer¬ 
sality. Although an analytical/approximate description 
of the postmerger dynamics is not available, the gauge- 
invariant Af,(j) curves contain, in analogy to the merger 
case, significant information about the system dynam¬ 
ics m- Specihcally, we interpret as being gener¬ 

ated by some Hamiltonian flow that continuously con¬ 
nects merger and postmerger. In terms of this Hamilto¬ 
nian evolution, the values provide 

initial conditions for the dynamics of the MNS/HMNS; 
it is then plausible to assume that the postmerger cor¬ 
relation follows from these initial conditions by conti¬ 
nuity. In order to assess this conjecture, we define the 
frequency given by the equation Mil = dEf,/dj, notice 
that and show that H is the relevant 

dynamical frequency for both inspiral-merger and post¬ 
merger. Recalling that the standard quadrupole formula 
predicts that a generic source with m = 2 geometry and 
rotating at frequency H emits GWs at a frequency 20, 
we plot the latter in Fig. and indeed observe that it 
corresponds to the main emission channel 0 J 22 during the 
whole evolution. In practice, the gauge-invariant O can 
be interpreted as the orbital frequency during the inspi¬ 
ral, and the angular frequency of the MNS/HMNS during 
postmerger. Furthermore, since merger remnants from 
larger K 2 binaries are less bound and have larger angular 
momentum support at formation, O(k^) (so / 2 ) must be 
a monotonically decreasing function of which is what 
one can observe in Fig. 

The frequency evolution is also expected to depend on 
angular momentum dissipation due to magnetic fields in¬ 
stabilities, e.g. [HHIS]) cooling and shear viscosity m- 
However, the available literature indicates these physical 
effects are negligible in first approximation, and we argue 
that they might result in frequency shifts A /2 < 6 / 2 - The 
stars rotation can instead play a relevant role via spin- 


orbit coupling effects: stars with dimensionless spin pa¬ 
rameters > 0.05 — 0.1 can give frequency shifts ^ 6/2 [39] . 

Outlook. — The result of this work, coupled with the 
modeling of the merger process given in mim, indi¬ 
cates the possibility to model the late-inspiral-merger- 
postmerger GW spectrum in a consistent way using as 
main parameter. In particular, an accurate late-inspiral- 
merger GW spectrum is given by a suitable frequency- 
domain representation, h{f) = A{f) exp [—i'I'(/)], of the 
waveform of m- The leading-order tidal contribution 
of such a spectrum reads = —39/4 k 2 with 

x{f) oc see [T] for at 2.5 post-Newtonian 

order. A simple template for the postmerger spectrum 
for binaries with M < Mthr is then given by a single¬ 
peak-model and our fit for / 2 . The precise construction 
of such complete spectrum will be subject of future work. 
As mentioned in our discussion, it will be particular im¬ 
portant to include spin effects, e.g. 13139]. 

The performance of the proposed model in a GW data- 
analysis context will be carefully evaluated in a separated 
study. In this respect, we suggest that an optimal strat¬ 
egy to constrain the EOS could be combining the late- 
inspiral measurement of type (I) with measurement of 
type (H). The inclusion of the postmerger model might 
lead to an improved estimate of for the same number 
of observed events |3H3. 
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